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We report a detection of the baryon acoustic oscillation (BAO) feature in the three-dimensional correlation function of the transmitted 
flux fraction in the Lya forest of high-redshift quasars. The study uses 48,640 quasars in the redshift range 2.1 < z < 3.5 from the 
Baryon Oscillation Spectroscopic Survey (BOSS) of the third generation of the Sloan Digital Sky Survey (SDSS-III). At a mean 
redshift z = 2.3, we measure the monopole and quadrupole components of the correlation function for separations in the range 
20h~'Mpc < r < 2001r 1 Mpc. A peak in the correlation function is seen at a separation equal to (1.01 ± 0.03) times the distance 
expected for the BAO peak within a concordance ACDM cosmology. This first detection of the BAO peak at high redshift, when the 
universe was strongly matter dominated, results in constraints on the angular diameter distance D A and the expansion rate H at z = 2.3 
that, by themselves, require the existence of dark energy. Combined with constraints derived from Cosmic Microwave Background 
(CMB) observations, this result implies H(z = 2.3) = (224 ± 8)km s~'Mpc~', indicating that the time derivative of the cosmological 
scale parameter a = H(z = 2.3)/(l + z) is significantly greater than that measured with BAO at z ~ 0.5. This demonstrates that the 
expansion was decelerating in the range 0.7 < z < 2.3, as expected from the matter domination during this epoch. Combined with 
measurements of Ho, one sees the pattern of deceleration followed by acceleration characteristic of a dark-energy dominated universe. 

Key words, cosmology, Lya forest, large scale structure, dark energy 
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1. Introduction 

Baryon acoustic oscillations (BAO) in the pre-recombination 
universe have striking effects on the anisotropies of the Cosmic 
Microwave Background (CMB) and o n the l arge scale struc- 
ture (LSS) of matter ( Weinberg et al.l d2012l) and references 
therein). The BAO effects were first see n in the series of 
peaks in the CMB angular power spectrum (de Bern ardis et all 
2000) Subsequently, t he BAO relic at redshif t z ~ 0.3 was 
seen (Eisenst ein et all l2005b ICole et all [2005) as a peak in 
the galaxy-galaxy correlation function at a co-moving distance 
corresponding to the sound horizon at recombination. Fo r the 
WMAP7 cosmological parameters (Koma tsu et all I201 if) , the 
expected comoving scale of the BAO peak is r s = 153 Mpc, 
with an uncertainty of m 1%. 

The BAO peak in the correlation function at a redshift z 
appears at an angular separation A9 = r s /(l + z)Da(z) and 
at a redshift separation Az = r s H(z)/c, where Da and H are 
the angular distances and expansion rates. Measurement of the 
peak position at any redshift thus constrains the combinations of 
cosmological parameters that determine r s H and r s /DA- While 
the possibility of measuring both combinations is beginning 
to be exploited (Chuang & Wang), 120121 : IXu et all l2012h . most 
present measurements have concentrated on the combination 
Dy = [(1 + z) 2 D 2 A czl 'H] 1 ^ ', which determines the peak position 
for an isotropic distribution of galaxy pairs and an isotropic clus- 
tering strength. The "BAO Hubble diagram", Dy/r s vs.z, now 
includes the Slo an Digital Sky Survey (SPSS ) measurement 
(lEisen stein et al., 2005) updated to the DR7 (Abaz aiian et all 
2009) sample and combined with 2dF data ( Percival et all 
20101) , the 6dF point at z = 0.1 (iBeutler et all 1201 ll) , the 
WiggleZ points at (0.4 < z < 0.8) (Bla ke et all 1201 1 al) . and 
a reanalysis of th e SDSS DR 7 sample that uses rec onstruction 
(lEisenstein et al. L 120071; iPadm anabha n et all 20091) to sharpen 
the p r ecision of th e BAO measurement (Padmana bhan et ail 
120 lit iMehta et all |2012|) . Re cently, the Baryon Oscillation 
Spe ctroscopic Survey (BO SS; iDawson et all 1201 2|) of SDSS- 
III (lEisenstein et all 201 ll) has added a precise measurement at 
z - 0.57 (lAnderson et alll2012l) . BOSS has also reported a mea- 
surement of Da(z = 0.55)/r s based on the BAO struc ture in the 
angular power spectrum of galaxies dSeo et alll2012l) . 

In this paper, we present an observation of the BAO peak at 
z ~ 2.3 found in the flux correlation function of the Lya- for- 
est of BOSS quasars. This is the first such observation at a red- 
shift where the expansion dynamics is matter-dominated, z > 
0.8. The possibilit y of such a measu rement was suggested by 
IMcDonald I d2003h" and White (2003) and first stud ed in detail 
bv lMcDonald & Eisenstein! (l2007l) Twhile the galaxy BAO mea- 
surements are most sensitive to Dy D 2 Ph~ 1 ^, the Lya- flux 
transmission is more sensitive to peculiar velocity gradient ef- 
fects, which enhance redshift distortions and shift our sensitivity 
to the expansion rate. As we shall show below, the most accu- 
rately measured combination from the Lyo- forest BAO peak is 
oc D^' 2 // -0 ' 8 , and the present BOSS data set allows us to de- 
termine its value to a precision of 3.5%. Combining this result 
with constraints from CMB observations allows us to deduce the 
value of H(z = 2.3) accurate to 4%. Comparing our results with 
measurements of Ho and of H(0.2 < z < 0.8) reveals the ex- 
pected sequence of deceleration and acceleration in models with 
dark energy. 

The last decade has seen increasing use of Lya absorption to 
investigate large sca le structure. The number of quasars in early 
studie s dCroft et all 119991: iMcDonald et all l2000t ICroft et all 



to determine the Lya absorption correlation along individual 
lines of sight. With the BOSS project the surface density of 
quasars is sufficient to probe the full three-dimensional distri- 
bution of neutral hydroge n. A study us i ng the first 10,000 BOSS 
quasars was presented bv lSlosar et al.l (1201 lh . This sample pro- 
vided clear evidence for the expected long-range correlations, 
including the redshift-space distortions due to the g ravitational 
grow th of structure. With the SDSS data release DR9 (lAhn et al. 
2012), we now have ~ 60,000 quasars at z ~ 2.3 (Pari s et al. 



2012), with a high enough surface density to observe the BAO 
peak. 

The use of Lya absorption to trace matter has certain inter- 
esting differences from the use of galaxies. Galaxy surveys pro- 
vide a catalog of positions in redshift space that correspond to 
points of high over-densities. On the other hand, the forest region 
of a quasar spectrum provides a complete mapping of the ab- 
sorption over a ~ 400 IT 'Mpc (comoving) range starting about 
lOOfr'Mpc in front of the quasar (so as to avoid the necessity 
of modeling the quasar's Lye emission line). To the extent that 
quasar lines of sight are random, a large collection of quasars 
can provide a nearly unbiased sample of points where the 
abso rption is measur ed. C osmological simula tions dCen et all 
" 1 994t iPetitiean et all 1 1 9951: IZhang et all 1 1 995b iHernquist et all 
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l2002HViel et alll2004H McDonald etalll2006h was enough only 



iMiralda-Escude et all 119961: iTheuns et all 119981) " indi- 
cate that most of the Lya absorption is due to cosmic fil- 
amentary structures with overdensities of order one to ten, 
much lower than the overdensities of virialized halos sam- 
pled by galaxies. These simulations have also indicated that, 
on large scales, the mean Lya absorption is a linear tracer of 
the mass overdensity (ICroft et all 1 19971 [19981 1 Weinberg et all 
1 19981: IMcDonald et allfeoOOt IMcDonald 1 120031) . implying a re- 
lation of the power spectrum of the measured absorption to that 
of the underlying mass fluctuations. Finally, the forest is observ- 
able in a redshift range inaccessible to current large galaxy sur- 
veys and where theoretical modeling is less dependent on non- 
linear effects in cosmological structure formation. These factors 
combine to make Lya absorption a promising tracer of mass that 
is complementary to galaxy tracers. 

The use of the Lya forest is handicapped by the fact that 
not all fluctuations of the transmitted flux fraction are due to 
fluctuations of the density of hydrogen. Because the neutral hy- 
drogen density is believed to be determined by photo-ionization 
equilibrium with the flux of U V photons from stars and quasars, 
variations of the UV flux may contribute to fl uctuations of the 
transmitted flux (IWorseck & Wisotzkil 120061) . In addition, hy- 
drogen can be self-shielded to ionizing photons, leading to much 
higher n eutral fractions and damping wings in strong absorption 
systems (iFont-R ibera et al., 2012b). Metals present in the inter- 
galactic medium pro vide additional ab sorption superimposed on 
the Lya absorption dPieri et all 120101) . A further complication 
results from the fact that the flux correlation function uses the 
fluctuations of the transmitted flux about its mean value; this re- 
quires an estimate of the product of the mean absorption (as a 
function of abs orber redshift) and t he unabsorbed flux for indi- 
vidual quasars dLe Goff et al. . l20llT) . 

All of these complications are important for a complete un- 
derstanding of the statistics of the Lya forest. Fortunately, they 
are not expected to produce a sharp peak-like feature in the cor- 
relation function, so the interpretation of a peak as due to BAO 
should give robust constraints on the cosmological parameters. 

The BOSS collaboration has performed three independent 
analyses to search for BAO in Lya forest. This paper presents 
two of them, both of which aim to analyze the forest using sim- 
ple procedures at the expense of some loss of sensitivity. The 
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third analysis, with the goal of an optimal measurement of the 
flux correlation function with a m ore complex method, will b e 
described in a separate publication dSlosar et aUlin preparation!) . 
Our methods are tested extensively on a set of detailed mock 
catalogs of the BOSS Lya forest dat a set. These mock catalogs , 
which use the method presented by iFont-Ribera et all J2012a), 
will be described in detail in a forthcoming public release paper 
(S. Bailey et al., in preparation). While our analyses implement 
their own sample selection criteria and continuum determination 
procedures, we have also applied our measurement to a "stan- 
dard" BOSS DR9 Lya forest sample that will be presented in a 
future paper (Lee et al., in preparation). 

This paper is organized as follows. Section [2] presents the 
BOSS quasar sample used in this analysis and the procedure 
used to produce the quasar spectra. Section[3]describes the anal- 
ysis to measure the correlation function. Section [4] derives the 
monopole and quadrupole components of the correlation func- 
tion and determines the significance of the peak observed in 
these functions at the BAO scale. The cosmological implica- 
tions of our detection of a BAO peak are discussed in Section 
[5] Finally, Appendix [A] provides a brief description of the mock 
spectra used to test our methodology and AppendixlBl shows the 
result of our BAO measurement applied to the BOSS Lya sam- 
ple of Lee et al. (in preparation). 

2. The BOSS quasar sample and data reduction 



The BOSS project dDawson et all 120121) of SDSS-III 
dEisenstein et all [201 ll) is obtaining the spectra of ~ 1.6 x 10 6 
luminous galaxies and ~ 150,000 quasars. The project uses 



upgraded versions of the SDSS spectrographs (Smeeetal 



2012) mounted on the Sloan 2.5-meter telescope (Gun n et al 



2006) at Apache Point, New Mexico. BOSS galaxy and quasar 
spectroscopic targets are selected using algori thms based 



1998 : lYork et all 



prima r ily on pho tometry from the SDSS camera (Gunn et al 
2000) in the ugriz bands (Fukugit a et al 
I2002h reduced and calibrated as de- 
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scribed by U toughton etalJ d2002h . iPier et all d2003h . and 
Padmanabhan et al. (2008). Targets are assigned to fibers ap- 
propriately positioned in the 3° diameter focal plane accordin g 
to a specially designed tiling algorithm dBlanton et all [2003). 
Fibers are fixed in place by a pierced metal plate drilled for 
each observed field and fed to one of two spectrographs. Each 
exposed plate generates 1,000 spectra covering wavelengths 
of 360 to 1000 nm with a r esolving power ranging from 1500 
to 3000 dSmee et all 120121) . A median of 631 of these fibers 
are assigned to galactic targets and 204 to quasar targets. The 
BOSS spectroscopic targets are observed in dark and gray time, 
while the bright-time is used by other SDSS-III surveys (see 
lEisensteinetal]l201ll) . 

The quasar spectroscopy targets are selecte d from photo- 
metri c data with a comb i nation of algorith ms (Richards et al. 
2009; Ye che et alJ 120091: iKirkpatrick et alJ l201lt iBovvetal. 
201 U iPal anque-De labrouille et all 1201 1[ for a summary, see 
Ross et all d2012l) I The algorithms use SDSS fluxes and, for 
SDSS Stripe 82, photometric variabilit y. When available , we 
also use data fro m non-optica l surve ys (Bow et al., 2012j): the 
GALE X survey dMartin et all l2005h in the UV; the UKIDSS 
survey (Lawrence et al,, 2007) in the NIR, and the FIRST sur- 
vey (Bec ker et all 1 1995b in the radio. 

The quasar spectroscopy targets are divided into two sam- 
ples "CORE" and "BONUS". The CORE sample consists of 20 
quasar targets per square degree selected from SDSS photome- 
try with a uniform algorithm, for which the selection efficiency 




Fig. 1. Hammer- Aitoff projection in equatorial coordinates of the 
BOSS DR9 footprint. The observations cover ~ 3000 deg 2 . 



for z > 2.1 qua sars is ~ 50%. Th e selection algorithm for the 
CORE sample (iBovv et all 1201 ll) was fixed at the end of the 
first year of the survey, thus making it useful for studies that 
require a uniform target selection across the sky. The BONUS 
sample was chosen from a combination of algorithms with the 
purpose of increasing the density on the sky of observed quasars 
beyond that of the CORE sample. The combined samples yield 
a mean density of identified quasars of 15 deg~ 2 with a maxi- 
mum of 20 deg -2 , mostly in zones where photometric variabil- 
ity, UV, and/or NIR data are available. The combined BONUS 
plus CORE sample can be used for Lya BAO studies, which re- 
quire the highest possible quasar density in a broad sky area but 
are insensitive to the uniformity of the quasar selection criteria 
because the structure being mapped is in the foreground of these 
quasar back-lights. 

The da t a pre sented here consist of the DR9 data release 
(lAhn et all 120121) covering ~ 3000 deg 2 of the sky shown in 
figureQ] These data cover about one-third of the ultimate BOSS 
footprint. 

The d ata were reduced w ith the SDSS-III pipeline as de- 
scribed in [Bolton et al.l d2012l) . Typically four exposures of 15 
minutes were co-added in pixels of wavelength width ~ 0.09 nm. 
Besides providing flux calibrated spectra, the pipeline provided 
preliminary object classifications (galaxy, quasar, star) and red- 
shift estimates. 

The spectra of all quasa r targe ts were visually inspected, 
as described in Pa ris et al. (2012), to correct for misiden- 
tifications or inaccurate redshift determinations and to flag 
broad absorption lines (BAL). Damped Lya troughs are visu- 
ally fl agged, b ut al s o ide ntified and characterized automatically 
dNoterdaeme et all 120121) . The visual inspection of DR9 con- 
firms 60,369 quasars with 2.1 < z q < 3.5. In order to simplify the 
analysis of the Lya forest, we discarded quasars with visually 
identified BALs and DLAs, leaving 48,640 quasars. 

For the measurement of the flux transmission, we use the 
rest-frame wavelength interval 



104.5 < A lf < 118.0nm . 



(1) 



The range is bracketed by the Lya and Ly/3 emission lines at 
121.6 and 102.5 nm. The limits are chosen conservatively to 
avoid problems of modeling the shapes of the two emission lines 
and to avoid quasar proximate absorbers. The absorber redshift, 
Z = A/Aiya - 1, is in the range 1.96 < z < 3.38. The lower 
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limit is set by the requirement that the observed wavelength be 
greater than 360 nm below which the system throughput is less 
than 10% its peak value. The upper limit comes from the max- 
imum quasar redshift of 3.5, beyond which the BOSS surface 
density of quasars is not sufficient to be useful. The distribution 
of absorber redshift is shown in figure[2](top panel). When given 
the weights used for the calculation of the correlation function 
(section [3~3l ). the absorbers have a mean redshift of (z) = 2.31. 

For the determination of the correlation function, we use 
"analysis pixels" that are the flux average over three adjacent 
pipeline pixels. Throughout the rest of this paper, "pixel" refers 
to analysis pixels unless otherwise stated. The effective width 
of these pixels is 210 km s _1 , i.e. an observed-wavelength width 
~ 0.27 nm ~ 2 fr'Mpc. The total sample of 48,640 quasars thus 
provides ~ 8 x 10 6 measurements of Lya absorption over a total 
volume of ~ 20fT 3 Gpc 3 . 

Figure [2] (bottom panel) shows the distribution of the signal- 
to-noise ratio for pixels averaged over the forest region. The rel- 
atively modest mean value of 5.17 reflects the exposure times 
necessary to acquire such a large number of spectra. 

In addition to the BOSS spectra, we analyzed 15 sets of 
mock spectra that were produced by the methods described in 
appendix [A] These spectra do not yet reproduce all of the char- 
acteristics of the BOSS sample, but they are nevertheless useful 
for a qualitative understanding of the shape of the measured cor- 
relation function. More importantly, they are useful for under- 
standing the detectability of a BAO-like peak and the precision 
of the measurement of its position. 

3. Measurement of the correlation function 

The flux correlation function can be determined through a sim- 
ple two-step process. In the first step, for each pixel in the forest 
region (equation [TJ of quasar q, the measured flux f q (A) at ob- 
served wavelength A is compared with the mean expected flux, 
Cq(A)F(z), thus defining the "delta field": 



5 q {A) 



fqW 

C„(A)F{z) 



1 



(2) 



Here, C q (A) is the unabsorbed flux (the so-called "continuum") 
and F(z) is the mean transmitted fraction at the HI absorber red- 
shift. The quantities A and z in equation [2] are not independent 
but related via z - A/A^ya - 1. 

Figure [3] shows an example of an estimation for C q {A) (blue 
line) and C q F (red line). Our two methods for estimating C q and 
F are described in sections [3~T1 and [3~2l 

In the second step, the correlation function is calculated as a 
weighted sum of products of the deltas: 



(3) 



ijeA 



ijeA 



where the wy are weights and each i or j indexes a measurement 
on a quasar q at wavelength A. The sum over (i, j) is understood 
to run over all pairs of pixels of all pairs of quasars within A 
defining a region in space of pixel separations, r, i-rj. The region 
A is generally defined by a range r m ; n < r < r max and ji m i a < /j < 
A<max with: 



(n-n) 



(4) 



where (r, - rf)\\ is the component along the line of sight. 
Separations in observational pixel coordinates (ra,dec,z) are 



wavelength (nm) 
360 400 440 480 520 




2.4 2.6 2.8 3.0 
redshift of absorbers 




5 10 

S/N per pixel 



20 



Fig. 2. Top: weighted distribution of absorber redshifts used in 
the calculation of the correlation function in the distance range 
80 Ir'Mpc < r < 120 h~'Mpc. Bottom: distribution of signal- 
to-noise ratio for analysis pixels (triplets of pipeline pixels) av- 
eraged over the forest region. 



transformed to (r, n) in units of h 'Mpc by using a ACDM fidu- 
cial cosmology with matter and vacuum densities of 



(Atf,Q A ) = (0.27,0.73) . 



(5) 



In the sum Q, we exclude pairs of pixels from only one 
quasar to avoid the correlated errors in 6, and Sj coming from 
the estimate of C q . Note that the weights in eq.[3|are set to zero 
for pixels flagged by the pipeline as having problems due, e.g., 
to sky emission lines or cosmic rays. 

A procedure for determining £ is defined by its method for 
estimating the expected flux C q F and by its choice of weights, 
Wij. The two methods described here use the same technique to 
calculate weights but have different approaches to estimate C q F. 
We will see that the two methods produce correlation functions 
that have no significant differences. However, the two indepen- 
dent codes were invaluable for consistency checks throughout 
the analysis. 

The two methods were "blind" to the extent that many of the 
procedures were defined during tests either with mock data or 
with the real data in which we masked the region of the peak 
in the correlation function. Among those aspects fixed in this 
way were the quasar sample, the continuum determination, the 
weighting, the extraction of the monopole and quadrupole cor- 
relation function and the determination of the peak significance 
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(section 0). This early freezing of procedures resulted in some 
that are suboptimal but which will be improved in future analy- 
ses. We note, however, that the procedures used to extract cos- 
mological information (section [5]l were decided on only after 
de-masking the data. 




440 460 



480 500 
X (nm) 



520 540 



Fig. 3. An example of a BOSS quasar spectrum of redshift 
3.239. The red and blue lines cover the forest region used here, 
104.5 < /Itf < 118.0. This region is sandwiched between the 
quasar's Ly/3 and Lya emission lines respectively at 435 and 
515 nm. The blue line is an estimate of the continuum (unab- 
sorbed flux) by method 2 and the red line is the estimate of the 
product of the continuum and the mean absorption by method 1 . 



3.1. Continuum fits, method 1 

Both methods for estimating the product C q F assume that C q is, 
to first approximation, proportional to a universal quasar spec- 
trum that is a function of rest-frame wavelength, Aa — A/( 1 + z q ) 
(for quasar redshift z q ), multiplied by a mean transmission frac- 
tion that slowly varies with absorber redshift. Following this as- 
sumption, the universal spectrum is found by stacking the ap- 
propriately normalized spectra of quasars in our sample, thus 
averaging out the fluctuating Lya absorption. The product C q F 
for individual quasars is then derived from the universal spec- 
trum by normalizing it to account for the quasar's mean forest 
flux and then modifying its slope to account for spectral-index 
diversity and/or photo-spectroscopic miscalibration. 

Method 1 estimates directly the product C q F in equation [2] 
An example is given by the red line in figure [5] The estimate is 
made by modeling each spectrum as 



(6) 



where a q is a normalization, b q a "deformation parameter", and 
(A) is the mean wavelength in the forest for the quasar q and 
/(/Irf , z) is the mean normalized flux obtained by stacking spectra 
in bins of width Az — 0. 1 : 



/(A rf ,z) = J> 9 / ? (W ? 128 /J/ 



'q ■ 



(7) 



Here z is the redshift of the absorption line at observed wave- 
length A (z - A/Aiya - 1), f q is the observed flux of quasar q 



at wavelength A and /_ 128 is the average of the flux of quasar 
q for 127.5 < Arf < 128.5 nm. The weight w q (A) is given by 
w~ 1 = l/[ivar(/l) • (/ ? 128 ) 2 ] + o~j. [ux LSS . The quantity ivar is the 
pipeline estimate of the inverse flux variance in the pixel corre- 
sponding to wavelength A. The quantity crL^ LSS is the contri- 
bution to the variance in the flux due to the LSS. We approxi- 
mate it by its value at the typical redshift of the survey, z ~ 2.3: 
^fiux lss ~ 0-035 (section [33T >. 

Figure [4] shows the resulting mean <5, as a function of ob- 
served wavelength. The mean fluctuates about zero with up to 
2% deviations with correlated features that include the H and K 
lines of singly ionised calcium (presumably originataing from 
some combination of solar neighborhood, interstellar medium 
and the Milky Way halo absorption) and features related to 
Balmer lines. These Balmer features are a by-product of imper- 
fect masking of Balmer absorption lines in F-star spectroscopic 
standards, which are used to produce calibration vectors (in the 
conversion of CCD counts to flux) for DR9 quasars. Therefore 
such Balmer artifacts are constant for all fibers in a plate fed 
to one of the two spectrographs and so they are approximately 
constant for every 'half-plate'. 

If unsubtracted, the artifacts in figure [4] would lead to spuri- 
ous correlations, especially between pairs of pixels with separa- 
tions that are purely transverse to the line of sight. We have made 
a global correction by subtracting the quantity (6)(A) in figure 0] 
(un-smoothed) from individual measurements of 6. This is justi- 
fied if the variance of the artifacts from half-plate-to-half-plate is 
sufficiently small, as half-plate-wide deviations from our global 
correction could, in principle add spurious correlations. 

We have investigated this variance both by measuring the 
Balmer artifacts in the calibration vectors themselves and by 
studying continuum regions of all available quasars in the DR9 
sample. Both studies yield no detection of excess variance aris- 
ing from these artifacts, but do provide upper limits. The study 
of the calibration vectors indicate that the square-root of the vari- 
ance is less than 20% of the mean Balmer artifact deviations and 
the study of quasar spectra indicate that the square-root of the 
variance is less than 100% of the mean Balmer artifacts (and 
less than 50% of the mean calcium line deviations). 

We then performed Monte Carlo simulations by adding a 
random sampling of our measured artifacts to our data to con- 
firm that our global correction is adequate. We found that there 
is no significant effect on the determination of the BAO peak po- 
sition, even if the variations are as large as that allowed in our 
tests. 



3.2. Continuum fits, method 2 

Method 1 would be especially appropriate if the fluxes had a 
Gaussian distribution about the mean absorbed flux, C q F. Since 
this is not the case, we have developed method 2 which explicitly 
uses the probability distribution function for the transmitted flux 
fraction F, P(F,z), where < F < 1. We use the P(F,z) that 
results from the log-normal model used to generate mock data 
(see appendix lAl. 

Using P(F,z), we can construct for each BOSS quasar the 
PDF of the flux in pixel i, f„ by assuming a continuum C q (A£) 
and convolving with the pixel noise, cr,-: 



PiifhCJAiUd* dFP(F,zdexp 



jr. 

Jo 



-(C q F-0 



2^ 



(8) 
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The continuum is assumed to be of the form 



CJA) = (a q + b q A)f(Ai) 



(9) 



where /(/irf) is the mean flux as determined by stacking spectra 
as follows: 



f(A d ) = J] w q (A d ) [f q (Aa )/f™] / J] 



(10) 



as in equation [7] except that here there is no redshift binning. 
The parameters a q and b q are then determined for each quasar 
by maximizing a likelihood given by 



L(C q ) = Y\Pi[fi,C q (Ai)] 



(11) 



Figure [3] shows the C q (A) estimated for a typical quasar (blue 
line). 

The last element necessary to use equation (0 is the mean 
transmitted flux fraction F(z). If P{F,z) derived from the mocks 
were the true distribution of the transmitted flux fraction, then 
F(z) could simply be computed from the average of this distri- 
bution. Since this is not precisely true, we determine F(z) from 
the data by requiring that the mean of the delta field vanish for all 
redshifts. The F(z) we obtain is shown in figure[5] The unphysi- 
cal wiggles in the derived F(z) are associated with the aforemen- 
tioned residuals in 6(A) for method 1 (figure[4]l. 

There is one inevitable effect of our two continuum estimat- 
ing procedures. The use of the forest data in fitting the contin- 
uum effectively forces each quasar to have a mean absorption 
near that of the mean for the entire quasar sample. This ap- 
proach introduces a spurious negative correlation between pix- 
els on a given quasar even when well separated in wavelength. 
This negative correlation has no direct effect on our measure- 
ment of the flux correlation function because we do not use pixel 
pairs from the same quasar. However, the physical correlation 
between absorption on neighboring quasars causes the unphysi- 
cal negative correlation for individual quasars to generate a neg- 
ative contribution to the correlation measured with quasar pairs. 
Fortunately, this distortion is a smooth function of scale so it can 
be expected to have little effect on the observability or position 
of the BAO peak. This expectation is confirmed by analysis of 
the mock spectra (section|5]l. 



redshift 
2.6 2.8 
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Fig. 4. The mean of 6(A) plotted as a function of observed wave- 
length (method 1). Systematic offsets from zero are seen at the 
2% level. The calcium lines (393.4,396.8 nm) is present. The 
features around the hydrogen lines Hy, 6 and e (434.1, 410.2, 
397.0 nm) are artifacts from the use of F-stars for the photocali- 
bration of the spectrometer. 
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Fig. 5. Mean transmitted flux fraction as a function of redshift 
obtained from the continuum fits with method 2. Data are shown 
in black, mock-000 in red and the input mean transmitted flux 
fraction in blue. 



3.3. Weights 

A discussion on the optim al use of weights for the L yg corre- 
lation function is found in iMcO uinn & White! d20TTb ." Here we 
simply choose the weights w/j so as to approximately minimize 
the relative error on £4 estimated with equation ©. In the ap- 
proximation of uncorrelated pixels, the variance of £4 is 



Var(&) = 



[Zi.jeAWij] 2 



(12) 



where the pixel variance, includes contributions from both 
observational noise and LSS. The signal-to-noise ratio is: 



(If 



<|a) 2 ^ »fi/"V;)' 
Var(f A ) ZijeA 



(13) 



Because of LSS growth and redshift evolution of the mean ab- 
sorption, the gij depend on redshift and we use the measured de- 
pendence of the Id correlation function (Mc Donald et all 2006) 



tij(z) = (i+ziy /2 (i+zjr%(zo) 



7/2, 



y ~ 3.8 



(14) 



Maximizing the signal-to-noise ratio with respect to w, ; this 
gives: 



a + Zi y /2 (i + zjY /2 



■2 £2 

jj 



(15) 



For this expression to be used, we require a way of estimating 
the We assume that it can be decomposed into a noise term 
and a LSS term ((Tlss)' 



2 pipeline,! 2 , \ 



Zi - Aj/A Lyo - 1 , 



(16) 
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Fig. 6. Top panel: Inverse total variance in bins of redshift 
as a function of the pipeline inverse variance. Bottom panel: 
Parameters of the fit: the LSS contribution <x L ss (left) and the 
pipeline correction factor r\ (right) as a function of redshift. The 
lines show fits to the data as explained in the text. 



Figure [7] shows an example of £(r, fj) for the r bin centered 
on 34 Ir'Mpc. The blue dots are the data and the red dots are 
the mean of the 15 mocks. The function falls from positive to 
negative values with increasing fi, as expected from redshift dis- 
tortions. The effect is enhanced by the deformation due to the 
continuum subtraction. 

Figure [8] presents £(r,/j) averaged over three bins in /j. A 
clear peak at the expected BAO position, r s = 105 Ir'Mpc, is 
present in the bin 0.8 < ju < 1.0 corresponding to separation 
vectors within 37° of the line-of-sight. The curves show the best 
fits for a ACDM correlation function, as described in section [5] 

The data were divided into various subsamples to search 
for systematic errors in £(r,fj). For example, searches were 
made for differences between the northern and southern Galactic 
cap regions and between higher and lower signal-to-noise ra- 
tio quasars. No significant differences were found in the overall 
shape and amplitude of the correlation function. We also verified 
that the BAO peak position does not change significantly when 
wavelength slices of Lya forest data are eliminated, in particular 
slices centered on the Balmer features in figure [4] The peak po- 
sition also does not change significantly if the subtraction of the 
mean 6 (figure |4]i is suppressed. 



CO 



JJ_P 
o 




where cr 2 pipdine ; = [ivar(C 9 F) 2 ] 1 is the pipeline estimate of the 
noise-variance of pixel i and 77 is a factor that corrects for a pos- 
sible misestimate of the variance by the pipeline. 

We then organize the data in bins of a 2 . .. . and redshift. 

° pipeline j 

In each such bin, we measure the variance of which serves as 
an estimator of for the bin in question. The two functions j](z) 
and o~^ ss (z) can then be determined by fitting equation ( [To*] ). 

These fits are shown in figure|6] The top panel shows that the 
measured inverse variance follows the inverse pipeline variance 
until saturating at the redshift-dependent LSS variance (shown 
on the bottom left panel). For z > 3, there are not enough pixel 
pairs to determine rj(z) and cr^ ss {z). In this high redshift range, 
we assumed r\ — 1 and extrapolated cr L ss(z) with a second- 
degree polynomial fit to the z < 3 data. 

3.4. fan) 

The procedure described above was used to determine £(r,fi) 
through equation [3] in r-bins of width 4 Ir'Mpc (centered at 
2,6,..., 198 Ir'Mpc) and in //-bins of width 0.02, (centered at 
0.01, 0.03, ... 0.99). The 50 x 50 r - fi bins have an average of 
6 x 10 6 terms in the sum (0 with an average nominal variance 
of £ for individual bins of (10 4 ) 2 as given by (eqnlTZll. 



Fig. 7. £(r,//) vs./u 2 for the bin centered on r — 34 Ir'Mpc. The 
red dots are the mean of the 15 mocks and the blue dots are the 
data. 



4. The Monopole and quadrupole 

The analysis of the correlation function was perfor med in the 
frame work of the standard multipole decomposition (Hamilton, 
1992). For each bin in r we fit a monopole (I = 0) and 
quadrupole {( = 2) to the angular dependence: 

f(r,/z) = J] ft(r)JM0 = Kb(r) -&(r)/2] + [3&(r)/2]// 2 (17) 

£=0,2 

where Pf is the ^-Legendre polynomial. We ignore the small 
and poorly determined £ = 4 term. This fit is performed using 
a simple x 1 minimization with the nominal variance (equation 
n~2T > and ignoring the correlations between bins. This approach 
makes the fit slightly sub-optimal. (Later, we will correctly take 
into account correlations between r-bins of the monopole and 
quadrupole.) We also exclude from this fit the portion // < 0.1 
to avoid residual biases due to correlated sky subtraction across 
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Fig. 8. £(r,//) averaged over 0.1 < < 0.5, 0.5 < /i < 0.8 and 
0.8 < n < 1. The curves give fits (section|5]) to the data imposing 
concordance ACDM cosmology. The BAO peak is most clearly 
present in the data for ji > 0.8. 



quasars; this has a negligible impact on the fits and, at any rate, 
there is little BAO signal at low fi. 

Figure [9] displays the monopole and quadrupole signals 
found by the two methods. The two methods are slightly offset 
from one another, but the peak structure is very similar. Figure 
[9] also shows the combination £o + 0.1^2 which, because of the 
small monopole-quadrupole anti-correlation (section 14. U . is a 
better-determined quantity. The peak structure seen in figure[8]is 
also present in these figures. 



Because of the continuum estimation procedure (sections 
13.11 and I3.21 >. we can expect that the monopole and quadrupole 
shown in figure|9]are deformed with respect to the true monopole 
and quadrupole. The most important difference is that the mea- 
sured monopole is negative for 60 fr'Mpc < r < 100 fr'Mpc 
while the true ACDM monopole remains positive for all r < 
130 fr'Mpc. The origin of the deformation in the continuum es- 
timate is demonstrated in appendix[A]where both the true and es- 
timated continuum can be used to derive the correlation function 
(figure lA.U . As expected, the deformation is a slowly varying 
function of r so neither the position of the BAO peak nor its am- 
plitude above the slowly varying part of the correlation function 
are significantly affected. 

4. 1 . Covariance of the monopole and quadrupole 

In order to determine the significance of the peak we must es- 
timate the covariance matrix of the monopole and quadrupole. 
If the fluctuations <5, in equation ([3} in different pixels were 
uncorrelated, the variance of £4 would simply be the weighted 
products of the fluctuation variances. This yields a result that is 
~ 30% smaller than the true correlation variance that we com- 
pute below. The reason is, of course, that the 5-pairs are corre- 
lated, either from LSS or from correlations induced by instru- 
mental effects or continuum subtraction; this effect reduces the 
effective number of pairs and introduces correlations between 
(r,fj) bins. 

Rather than determine the full covariance matrix for £(r, //), 
we determined directly the covariance matrix for %o(r) and £2(7") 
by standard techniques of dividing the full quasar sample into 
subsamples according to position on the sky. In particular we 
used the sub-sampling technique described below . We also tried 
a bootstrap technique de.g. Efron & Gongl Il983h consisting of 
substituting the entire set of N subdivisions of the data by N 
of these subdivisions chosen at random (with replacement) to 
obtain a "bootstrap" sample. The covariances are then measured 
from the ensemble of bootstrap samples. Both techniques give 
consistent results. 

The adopted covariance matrix for the monopole and 
quadrupole uses the sub-sampling technique. We divide the data 
into angular sectors and calculate a correlation function in each 
sector. Pairs of pixels belonging to different sectors contribute 
only to the sector of the pixel with lower right ascension. We 
investigated two different divisions of the sky data: defining 800 
(contiguous but disjoint) sectors of similar solid angle, and tak- 
ing the plates as defining the sectors (this latter version does not 
lead to disjoint sectors). The two ways of dividing the data lead 
to similar covariance matrices. 

Each sector s in each division of the data provides a mea- 
surement of ^ s (r,fi) that can be used to derive a monopole and 
quadrupole, & s (r), (t = 0, 2). The covariance of the whole BOSS 
sample can then be estimated from the weighted and rescaled co- 
variances for each sector: 



= ( ^W s (r)W s (r') [hs{r)l e , s {r') - l t {r%{r')] ) 



(18) 



The average denoted by ( ) is the simple average over sec- 
tors, while £ f (r) denotes the correlation function measured for 
the whole BOSS sample. The W s (r) are the summed pixel-pair 
weights for the radial bin r for the sector s and W(r) is the same 
sum for the whole BOSS sample. 

The most important terms in the covariance matrix are the 
r — r' terms, i.e. the monopole and quadrupole variances. They 
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Fig. 9. Monopole (upper panel) and quadrupole (middle panel) 
correlation functions found by method 1 (red) and method 2 
(black). The bottom panel shows the combination £o + 0.1^2 
found by method 1 (red) and method 2 (black). 



are shown in figure [10] as a function of r. In the figure, they are 
multiplied by the number N of pixel pairs in the r-bin. The prod- 
uct is nearly independent of r, as expected for a variance nearly 
equal to the pixel variance divided by N. For the monopole, the 
variances are only about 30% higher than what one would cal- 
culate naively assuming uncorrelated pixels and equation (fT2l) . 
Figure [10] also displays the monopole-quadrupole covariance 
times number of pairs, which also is nearly independent of r. 

FigureQT]displays the monopole-monopole and quadrupole- 
quadrupole covariances. Nearest-neighbor covariances are of or- 



der 20%. Figure [TT] also shows monopole-quadrupole covari- 
ance. 

We used the 15 sets of mock spectra to test our method for 
calculating the covariance matrix. From the 15 measurements 
of ^c(r) one can calculate the average values of ^(r)^(r') and 
compare them with those expected from the covariance ma- 
trix. Figures Q~2] shows this comparison for the monopole and 
quadrupole variance, the monopole and quadrupole covariances 
between neighboring r-bins and the monopole-quadrupole co- 
variance. The agreement is satisfactory. 

4.2. Detection significance of the BAO peak 

In this section, we estimate the significance of our detection of a 
BAO peak at 105 h _1 Mpc. At the statistical power of the present 
data, it is clear that the peak significance will depend to some 
extent on how we treat the so-called "broadband" correlation 
function on which the peak is superimposed. In particular, the 
significance will depend strongly on the r-range over which the 
correlation functions are fitted. To the extent that the BAO peak 
is known to be present in the matter correlation function and that 
the Lya absorption is known to trace matter, the actual signif- 
icance is of limited interest for cosmology. Of greater interest 
is the uncertainty in the derived cosmological parameter con- 
straints (section [5]l which will be non-linear reflections of the 
peak significance derived here. 

A detection of the BAO peak requires comparing the quality 
of a fit with no peak (the null hypothesis) to that of a fit with a 
peak. Typically, this exercise would be performed by choosing 
a test statistic, such as the^ 2 , computing the distribution of this 
quality indicator from a large number of peak-less simulations 
and looking at the consistency of the data with this distribution. 
Since our mock data sets are quite computationally expensive 
and only a handful are available, we chose a different approach. 

Our detection approach uses the following expression to fit 
the observed monopole and quadrupole. 



B4V) + C4 ea V)+A^'(r) 



(19) 



where ^ B is a broadband term to describe the LSS correlation 

function in the absence of a peak, ^ eak is a peak term, and ^ lst is 
a "distortion" term used to model the effects of continuum sub- 
traction. The broadband term is derived from the fiducial ACDM 
cosmology defined by the parameters in equation (IA.U . It is ob- 
tained by fitting the shape of the fiducial correlation function 
with an 8 -node spline function masking the region of the peak 
(80 h'Mpc < r < 120 h"'Mpc). The peak term is the difference 
between the theoretical correlation function and the broadband 
term. Finally, the distortion term is calculated from simulations, 
as the difference in the monopole or quadrupole measured using 
the true continuum and that measured from fitting the continuum 
as described in appendix[A] The three components are shown in 
figure [T3l 

Expression ( TT9b contains three parameters each for the 
monopole and quadrupole (so six in total). We have performed 
fits leaving all six parameters free and fits where we fix the ra- 
tio C2/C0 to be equal to its nominal value used to generate our 
mock spectra (the value given by assuming a "redshift distor- 
tion parameter"/? = 1.4, see appendix [A}. We define the test 
statistic as the^- 2 difference between fitting equation ( fT9l simul- 
taneously to monopole and quadrupole by fixing Co to zero (a 
"peak-less" four or five-parameter fit) and fitting for Co (a five 
or six-parameter fit). In our detection fits we do not fit for the 
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Fig. 10. The r-dependence of the product of the monopole (top) 
and quadrupole (middle) variances and the number of pairs in 
the r-bin. The bottom panel shows the product for the monopole- 
quadrupole covariance (r = r' elements). The dotted lines show 
the means for r > 20fT I Mpc. 



BAO position but fix it to the theoretical prediction. The distri- 
bution for this test statistic ("Ay^ et ") under the null hypothesis is 
&X 2 distribution with one degree of freedom. The significance is 
then given by cr = (A^ et ) 1/2 . 

Figure [14] shows the fits to monopole (top panel) and 
quadrupole (bottom panel) and the corresponding fits with and 
without peaks and fixing ft = 1.4. For method 2, we obtain 
X 2 /DOF = 93.7/85 (111.8/86) with (without) a peak, giving 
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Fig. 11. The monopole and quadrupole covariance matrix. 
The monopole-monopole and quadrupole-quadrupole elements 
are normalized to the variance: Cy/ JCyCjj. The monopole- 
quadrupole elements are normalized to the mean of the 
quadrupole and monopole variances. The first off-diagonal el- 
ements of the monopole-monopole and quadrupole-quadrupole 
elements are ~ 20% of the diagonal elements. The diagonal ele- 
ments of the quadrupole-monopole covariance are 0.2 times 

the geometric mean of the monopole and quadrupole variances. 



Ax 2 = 18.1 for a detection significance of 4.2<x. For method 
1, we obtain x 2 /DOF = 93.2/85 (102.2/86) for a significance 
of 3. Oct. Allowing ft to be a free parameter gives essentially the 
same detection significances. 
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Fig. 12. Verification of the off-diagonal elements of the covari- 
ance matrix with the 15 sets of mock spectra. The black lines 
show correlations derived from the dispersion of the 15 mea- 
surements and the red lines show the correlations expected from 
the covariance matrix calculated by sub-sampling. The top and 
middle panels show the correlation between neighbouring bins 
for monopole and quadrupole respectively. The bottom panel the 
correlation between monopole and quadrupole measured at the 
same distance bin. 



The detection significance of ~ 4<x is typical of that which 
we found in the 15 sets of mock spectra. For the mocks, the 
significances ranged from to 6cr with a mean of 3.5cr. 

Our significance depends strongly on the fitting range. For 
a lower boundary of the range of r m ; n = 20,40,60 fr'Mpc 



Fig. 13. The fitting functions used for the determination of the 
peak detection significance: r 2 £ hb (r) (blue), r \^ eak {r) + £ bb (r)] 
(black) and r 2 ^ d,st (r) (red) for the monopole (solid lines) and 
quadrupole (dashed lines). 



we obtain a significance of <x = 4.2, 3.2 and 2.3, respectively 
(method 2). The reason for this result is illustrated in figure [T31 
where the results of the fits with and without peaks are com- 
pared to data for different values of r m j n = 20,40,60 fr'Mpc. 
Reducing the fitting range poses less stringent constraints on 
the distortion and broadband terms, thus allowing some of the 
peak to be attributed to the broadband. In particular, the statis- 
tically insignificant bump in the quadrupole at ~ 65 fr'Mpc 
causes the fitted broadband to increase as r m j n is increased to 
60 fr'Mpc, decreasing the amplitude of the BAO peak. For the 
monopole, the r m ; n = 60 fr'Mpc fit predicts a positive slope 
for £o(r) that decreases the amplitude of the peak but predicts a 
%o(r < 50 fr'Mpc) to be much less than what is measured. 

5. Cosmology with the BAO peak 

The observed position of the BAO peak in t;(r,/j.) is determined 
by two sets of cosmological parameters: the "true" parameters 
and the "fiducial" parameters. Nature uses the true cosmology to 
create correlations at the true sound horizon, r s . The true cosmol- 
ogy transforms physical separations between Lya absorbers into 
angles on the sky and redshift differences: Obao = r s / Da(z){\+z) 
and Azbao = f s H(z)/c. We, on the other hand, use a "fiducial" 
cosmology (defined by equation I A. U to transform angular and 
redshift differences to local distances at the redshift in question 
to reconstruct £(r,/j). If the fiducial cosmology is the true cos- 
mology, the reconstructed peaks will be at the calculated fiducial 
sound horizon, r s j. Limits on the difference between the fiducial 
and reconstructed peak position can be used to constrain the dif- 
ferences between the fiducial and true cosmological models. 

5.1. The peak position 

The use of incorrect fiducial Da, H and r s leads to shifts in the 
BAO peak position in the transverse and radial directions by the 
multiplicative factors a, and a r : 



D A (z)/r s 
D A,f/r s ,f 

Hf(z)r sJ 
H(z)r s 



(20) 



(21) 
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Fig. 14. Monopole and quadrupole fits with a BAO peak (red 
line) and without a BAO peak (blue line, Method 2). The fitting 
range is r m \ n < r < 200 fr'Mpc with r ra j n = 20 fr'Mpc. 



where th e subs cript / refers to the fiducial model. Following 
IXu et al.l (1201 2l) . we will use a fitting function, ^e(r), for the 
monopole and quadrupole that follows the expected peak posi- 
tion as a function of (a r , a r ): 

Ur) = Ur, at, a r , b,/3) + A e (r) { = 0,2. (22) 

Here, the two functions ff, derived from the power spectrum 
given in appendix [A] describe the underlying mass correlation 
function, the linear bias b and redshift distortion parameter 
and the movement of the BAO peak for a,, ay + 1. The func- 
tions At take into account distortions, as described below. 

For a t = a r = a, JO , there is a simple isotropic scaling 
of the coordinates by a, JO and |(r) is given by ff(r, a, S£) ) = 
ft(b,p)^(j(ai SO r), where £,tj are tne fiducial monopole and 
quadrupole and the normalizations ft are th e functions of th e 
bias and reds hift-distortion p arameter given by Hami lton! d 19921) . 
For a, + a r , IXu et al.l (1201 2l) found an approximate formula for 
|(r) that was good in the limit \a, - a r \ <s 1. We take the more 
direct route of numerically expanding £/(a,r f , a, r, ) in Legendre 
polynomials, PfQi), to directly calculate the |^(r, a t , a r ). 

The functions A[(r) describe broadband distortions due to 
continuum subtraction and the fact that the broadband correla- 
tion function is not expected to change in the same way as the 
BAO peak position when one deviates from the fiducial model. 



0.5 




Fig. 15. Same as figure [14] except for different fitting ranges: 
blue, green and red curves are for fits with r m ; n = 20 fr'Mpc, 
40 Ir'Mpc and 60 h _1 Mpc respectively (Method 2). The solid 
lines for fits without a BAO peak and the dashed lines with a 
peak. 

They correspond to the term A^ lst in equation ( fT9l . We have 
used two forms to represent A[(r): 

A f\ r) = q + h +Ce (23) 

and 

A f {r) = a S 2 + b -l+ Ci+ lL. (24 ) 

The observed monopole and quadrupole can then be fit to (1221 ) 
with free parameters a,, a r , bias, ft, and the nuisance parameters 
{a t , bt, c e and d e ). 

We first fixed {a,, a r ) = (1, 1) to determine if we find reason- 
able values of (b,/3). These two parameters are highly degenerate 
since both the quadrupole and monopole have amplitudes that 
are proportional to b 2 times polynomials in p. A well-determined 
combination is b(l + /3), for which we find a value 0.38 + 0.07; 
this is in agreemen t with b(l + J3) = 0.336 + 0.012 found at 
r ~ 40 fr'Mpc by Slosar et aT| (1201 ll) . The larger error of our 
fit reflects the substantial freedom we have introduced with our 
distortion function. 

We next freed all parameters to constrain (a,, ay). The con- 
tours for the two methods and two broadbands are shown in fig- 
ure [16] and the^- 2 for the fiducial and best-fit models are given in 
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Fig. 16. The contours for (D^/r^ r s H) obtained by fitting the 
monopole and quadrupole to (1221 . The broadband distortions are 
eqn. d23l dashed lines) or (l24l solid lines). The blue lines are for 
method 1 and the red lines for method 2. All contours are for 
Ax 2 = 4 except for the interior solid red contour which is for 
A* 2 = 1. 



table Q] The broadband term in equation (l24l) fits the data better 
than that in equation d23l both for the fiducial parameters and for 
the best fit. For broadband in equation (1241) . the x 2 f° r the fidu- 
cial model is acceptable for both methods: x 2 /DOF = 85.0/80 
for method 1 and x 2 /DOF = 71.5/80 for method 2. 

The contours in the figure are elongated along the direction 
for which the BAO peak position stays approximately fixed at 
large fi (near the radial direction, where the observations are 
most sensitive). The best constrained combination of Da and H 
of the form (D^H^ 1 /r s ) turns out to have f ~ 0.2. This low 
value of £ reflects the fact that we are mostly sensitive to the 
BAO peak in the radial direction. At the one standard-deviation 
level, the precision on this combination is about 4%. However, 
even this combination is sensitive to the tails in the contours. 
A more robust indicator of the statistical accuracy of the peak- 
position determination comes from fits imposing a, = a r = a, JO , 
as has gene rally been done in pre vious B AO studies wi th the 
exception of IChuang & Wangl d2012l) and lXu etalJ d2012l) . This 
constraint does not correspond to any particular class of cosmo- 
logical models. It does however eliminate the tails in the con- 
tours in a way that is similar to the imposition of outside data 
sets. The two methods and broadbands give consistent results, 
as seen in table Q] 

We used the sets of mock spectra to search for biases in 
our measurement of ff, so . The mean value reconstructed for this 
quantity on individual mocks is 1.002 ± 0.007, suggesting that 
there are no significant biases in the determination of the BAO 
scale. Figure [17] shows the values and errors for the individual 
mocks along with that for the data. Both the measured value and 
its uncertainty for the data is typical of that found for individual 
sets of mock spectra. 

5.2. Constraints on cosmological models 

Our constraints on (D^/r,, Hr„) can be used to constrain the cos- 
mological parameters. In a ACDM cosmology, apart from the 
pre-factors of Ho that cancel, Da/t s and Hr s evaluated at z — 2.3 
depend primarily on Q« through r s and on Qa which, with Qm, 
determines Da and H. The sound horizon also depends on Hq 



Fig. 17. The measurements of a, TO (= a, = a r ) for the 15 sets of 
mock spectra and for the data (realization=-l). The large errors 
for realization 5 and 8 are due to the very low significance of the 
BAO peak found on these two sets. 



Model: Open ACDM 




Fig. 18. Constraints on the matter and dark-energy density 
parameters (Qm, Qa) assuming a dark-energy pressure-density 
ratio w = — 1. The blue regions are the one and two stan- 
dard deviation constraints derived from our contours in figure 
[T6l (me thod 2, broadband l24l > combined with a measurement 
of Hq dRiess et all 1201 ll) . Also shown are one and two stan- 
dard deviation contours from lower redshift mea surements of 
Dy/r s (also combined with Hp) at z = 0.11 r 6dF: [Beutler et all 
(1201 lh l. z = 0.35 [LRG: | Percival et all (1201 Oh l and z = 
0.57 [CMASS: Anders on et alJ d2012l) l. All constraints use a 
WMAP7 (Komatsu et all 12011b prior on the baryon-to-photon 
ratio j] but do not otherwise incorporate CMB results. 

(required to derive Q r from Tqmb), on the effective number of 
neutrino species A^ v (required to derive the radiation density from 
the photon density), and on the baryon-to-photon number ratio, 
77 (required for the speed of sound). 

Figure [18] shows the ACDM constraints on (Qm, Qa) de- 
rived from the contours in figure [T6l combined with the most 
recent measurement of Hq dRiess et all 1201 ll) . We use the con- 
tours for method 2 and the broadband of equation [24] which 
gives better fits to the data than the other method and broadband. 
The contours also a ssume N v = 3 and the WMAP7 value of 77 
dKomatsu et all 1201 ll) . A lso shown are const r aints from BAO 
measurements of Dyjr, dPercival et all 1201 Ot lAnderson et all 
l2012HBeutler et all 1201 ll) . 
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Table 1. Results with the the two methods and two broadbands (equations l23l and l24l . Columns 2 and 3 give the^ 2 for the fiducial 
model and for the model with the minimum^- 2 . Column 4 gives the best fit for a (JO with the constraint (a, = ay = o-, J0 ). Column 5 
gives Hr J\Hr s ] gd with the 2 <x limits in parentheses. Column 6 gives the Hr s /[Hr s ] deduced by combining our data with that of 
WMAP7 (Komatsu e tafll201ll) (see section EDI 



method & 
broadband 


X)JD0F 


XIJDOF 


"(SO 


Hr s /[Hr s ] fid 


Hr s /[Hr s ] f u 
(with WMAP7) 


Method 1 l[24j 


85.0/80 


84.6/78 


1.035 ±0.035 


0.876 ± 0.049 Cg;lff) 


0.983 ± 0.035 


Method 2 (24) 


71.5/80 


71.4/78 


1.010 ±0.025 


0.954 ± 0.077 


1.000 ±0.036 


Method 1 ([23} 


104.3/82 


99.9/80 


1.027 ±0.031 


0.869 ± 0.044 Cg;J^) 


0.988 ± 0.034 


Method 2 i[23j 


88.4/82 


87.7/80 


1.004 ±0.024 


0.994 ±0.111 


1.006 ±0.032 



Model: Flat wCDM 



is -1.0 




Fig. 19. As in figure[T8]with constraints on the matter density pa- 
rameter, Qm, and dark-energy pressure-density ratio w assuming 
Qm + Qa — 1 . 



The Lya contours are nicely orthogonal to the lower redshift 
Dv/r s measurements, reinforcing the requirement of dark en- 
ergy from BAO data. In fact, our measurement is the only BAO 
measurement that by itself requires dark energy: Qa > 0.5. This 
is because at z = 2.3 the universe is strongly matter dominated 
and the VQm factor in H partially cancels the 1/ VQm in r s , 
enhancing the importance of the Qa dependence of H. 

Figure [19] shows the constraints on (Om,w; where w is the 
dark-energy pressure-density ratio) assuming a flat universe: 
Qa = 0. Our result is the only BAO measurement that by itself 
requires negative w. Our limit w < -0.6 requires matter domina- 
tion at z — 2.3. 



Pdeil = 2.3) 
pjz = 2.3) 



< 0.3 



Qa/Qm 
0.73/0.27 



(25) 



5.3. Constraints on H(z) 



The contours in figure[T6]give the measurements of Hr s given in 
table[T] A measurement of the expansion rate deep in the matter- 
dominated epoch can be used to demonstrate the deceleration 
of the expansion at that time. Unfortunately, our data is not yet 
precise enough to do this. To make a more precise measurement 
of H(z = 2.3), we must add further constraints to eliminate the 
long tails in figure [16] These tails correspond to models where 
l/H(z = 2.3) is increased (resp. decreased) with respect to the 
fiducial value while Da(z - 2.3) is decreased (resp. increased). 
For flat models, this requires the (l/H) (averaged up to z — 2.3 



to change in the opposite direction to that for l/H(z = 2.3), 
requiring an H(z) very far from fiducial. It is also possible to 
construct such models by introducing curvature. 

Because of the importance of curvature, the tails a re elim- 
inated once WMAP7 constraints ( Komats u et al.L 1201 ll) are in- 
cluded. This is done in figure [20] within the framework of non- 
flat models where the dark-energy pressure-density ratio, w (z), is 
determined by two parameters, wq and w a : w(z) = wq + w a z/(l + 
z). As expected, the WMAP7 results in this framework constrain 
Da and l/H to migrate in roughly the same direction as one 
moves away from the fiducial model. Combining WMAP7 con- 
straints with ours gives the values of H(z = 2.3)r 5 given in the 
last column of table [T] For what follows, we adopt the mean of 
methods 1 and 2 that use the more flexible broadband of equation 
C3J: 



H(z = 2.3)r s 
[H(z=23)r s ] fid 



0.992 ± 0.035 



(26) 



The precision on H is now sufficient to study the redshift evolu- 
tion of H(z). 

The fiducial model has r s = 152.76 Mpc and H(z - 2.3) = 
3.23//o, Hq = 70km s^'Mpc -1 . These results produce 



H(z = 2.3)r s 4 , 

— — — = (1.036 + 0.036) x 10 4 km s" 1 

1 +z 

or equivalently 



H(z = 2.3) 
1 + z 



(67.8 + 2.4)kms" 1 Mpc^ 1 



152.76 



(27) 



(28) 



This number can be compared with the measurements at 
lower redshift summarized in table [2] and shown in figure [2T] 
Other than measurements of Hq, the measurements of H(z) can 
be divided into two classes: those (like ours) that use r s as the 
standard of length and those that use c/Hq as the standard of 
length. Not shown are the less precise measurements tha t use 
stellar-population ev olution as a clock (IStern et all I20T0I) and 
the measurement of iRiess et al ] (I2007h that gave a model de- 
pendent estimate of H(z ~ 1 .2) by using the type la supernova 
Hubble diagram. 

The comparison with our measurement is simplest with mea- 
surements that use r s as the standard of length and therefore mea- 
sure H(z)r s (as is done here). There are four gro ups that have 
made this measurement. Chuang & Wand d2012l) and IXu et al.l 
(1201 2h studied the SDSS DR7 LRG sample and decomposed 
the BAO peak into radial and angular co mponents, thus ex- 
tracting directly H r s and DA/r s . Two groups (Bl ake et "all 120121 : 
iReid et all l2012h found H(z) by using the angle-averaged peak 
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Table 2. Measurements of H(z)/(l + z). The first seven mea- 
surements use r s = 1 52.75 tr'Mpc as th e standard of length. 
The measurements of iBlake et al.l d201 lbl) use supernova data 
and therefore measure H(z) relative to Hq. We quote the results 
they obtain without assuming a flat universe. 



2.3 

2.3 

0.35 

0.35 

0.57 

0.44 

0.60 

0.73 

0.2 
0.4 
0.6 
0.8 





H(z)/(l + z) 
kms _1 Mpc _1 



66.5: 
67.8: 
60.8: 
62.5: 
58.8: 
57.4: 
54.9: 
56.2: 

H X (1.11 : 

H x (0.83 : 

H X (0.81 : 

H x (0.83 : 



7.4 
2.4 
3.6 
5.2 
2.9 
5.4 
3.8 
4.0 



0.17) 
0.13) 

: 0.08) 

0.10) 



73.8 ±2.5 



method 



BAO 
BAO + WMAP7 

BAO 

BAO 
BAO + AP 
BAO + AP 



AP + SN 



reference 



this work 
this work 
Chuan g & Wang (2012) 
Xu et al. (2012) 
Reid etal. (2012) 
Blake ct al. (2012) 



1.2 
1.0 

0.8 

0.6 
0.7 









































WMAP 


wn6\ AD 


Lvinan-a 
^Lyman-a 









0.8 



0.9 1.0 

(Hxr s )/(Hxr s ) fid 



1.1 



1.2 



Fig. 20. Constraints on (DA/r s ,r s H) z= 2.3 within the frame- 
Blake et aL QOllh) work of OwOwaCDM models. The green contours are our lcr 
and 2cr constraints using method 2 and broadband (l24l . The 
gray contours a re the lcr and 2<x constraints from WMAP7 
( Ko matsu et al.L 1201 ll) . The red contours show the combined 

, ™,,s constraints. 
Riess etal. (2011) 



position to derive Dy(z)/r s = ((1 + z) \D A czH /r s . They then 
studied the Alcock-Paczynski effect on the broadband galaxy 
correlation function to determine Da(z)H(z). Combining the two 
measurements yielded H(z)r s . 

To demonstrate deceleration in the redshift range 0.5 < z < 
2.3, we fit the BAO-based values of H(z) in table|2]to the ACDM 
form H(z) = H (£l A + Q M (1 + z) 3 + (1 - Cl A - Q M )(1 + z) 2 ) 1 ' 2 . 
Marginalizing over Oa and Hq we find 



[H(Z)/(1+Z)] Z= 2.3 
[H(z)/(l+Z)] z=0 . 5 



= 1.17 ±0.05 , 



(29) 



clearly indicating deceleration between z = 2.3 and z = 0.5. 
This measurement is in good agreement with the fiducial value 
of 1.146. We emphasize that this result is independent of r s . 

To map the expansion rate for < z < 2.3, we must adopt the 
fiducial value of r s and compare the resulting H(z) with Hq and 
with ot her BAO-free mea surements. Besides the Hq measure- 
ment of Ri ess et a D (l2"()Tlt) . we use the WiggleZ analysis com- 
bining their Alcock-Paczynski data with distant super nova data 
from the Union-2 compilation dAmanullah et al.L l2010t) . The su- 
pernova analysis does not use the poorly known mean SNIa lu- 
minosity, so the SNIa Hubble diagram gives the luminosity dis- 
tance in units of Hr 1 , Dl(z)Hq. Combining this result with the 
Alcock-Paczynski measurement of Da(z)H(z) yields H(z)/Hq. 
The values are given in table|2] 

We fit all the data in table|2]to the ACDM form of H(z). This 
yields an estimate of the redshift of minimum H(z)/(1 + z) 




Zd-a = 0.82 ± 0.08 



(30) 



which compares well with the fiducial value: Zd-a - 
(2Q A /Q M ) 1/3 - 1 = 0.755 



6. Conclusions 

In this paper, we have presented the first observation of the BAO 
peak using the Lya forest. It represents both the first BAO de- 
tection deep in the matter dominated epoch and the first to use a 



Fig. 21. Measurements of H(z)/(1 + z) vs z demonstrating 
the acceleration of the expansion for z < 0.8 and deceleration 
for z > 0.8. The BAO-based measurements are our point at 
z = 2.3 (red square), the SDSS DR7 BAO me asurement at z = 
0.35 llXu et all (120121) : filled black square] and llChuang & Wanpl 
(120121) : open black square ! , the BOS S BAO - Alcock-Paczynski 
measurement at z = 0.57 llReid et alj d2012l) . cyan square!, and 
the W iggleZ BAO- Alcock-Paczynski measurement flBlake et al.l 
(20 1 2j), green squa re!. The green circles are from WiggleZ 
(Blak e et aLl 1201 lbh Alcock-Paczynski data combined with su- 
pernova data yielding H(z)/Hq (without the flatness assumption) 
plotted here assuming Hq — 70kms~ 1 Mpc~ 1 . The blue circle 
is the Hq measurement of iRiess et al.1 d201 lb . The line is the 
ACDM prediction for (h, Q. M , Q A ) = (0.7, 0.27, 0, 73). 



tracer of mass that is not galactic. The results are consistent with 
concordance ACDM, and require, by themselves, the existence 
of dark energy. Combined with CMB constraints, we deduce the 
expansion rate at z = 2.3 and demonstrate directly the sequence 
of deceleration and acceleration expected in dark-energy domi- 
nated cosmologies. 
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BOSS continues to acquire data and will eventually produce 
a quasar sample three times larger than DR9. We can thus ex- 
pect improved precision in our measurements of distances and 
expansion rates, leading to improved constraints on cosmologi- 
cal parameters. The Lya forest may well be the most practical 
method for obtaining precise Da{z) and H(z) measurements at 
z > 2, thanks to the large number of independent density mea- 
surements per quasar. It is reassuring that the first sample large 
enough to yield a detection of BAO produces a signal in good 
agreement with expectations. In the context of BAO dark energy 
constraints, high redshift measurements are especially valuable 
for breaking the degeneracy between curvature and the equation 
of state history More generally, however, by probing an epoch 
largely inaccessible to other methods, BAO in the Lya forest 
have the potential to reveal surprises, which could provide criti- 
cal insights into the origin of cosmic acceleration. 
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Appendix A: Mock quasar spectra 

We have produced mock spectra in order to tune the analysis 
procedure and to study statistical uncertainties and systematic 
effects in the measured correlation function. The basic process 
for gen erating these mocks is described in iFont-Ribera et alj 
(I2012ah . 

The mock quasars were generated at the angular positions 
and redshifts of the BOSS quasars. The unabsorbed spectra 
(continua) of the quasars were generated using the Principal 
Comp onent Analysis (PCA) eigenspectra of Suzuk i et alj 
(2005), with amplitudes for each eigenspectrum randomly drawn 
from Gaussian distributions wi th sigma equa l to the correspond- 



ing eigenvalues as published in Suzuki ( 2006) table 1 . A detailed 



description will be provided by Bail ey et al. 



in preparationl) . ac- 



companying a public release of the mock catalogs. 

We generated the field of transmitted flux fraction, F, that 
have a ACDM power spectrum with the fiducial parameters 



(fl«, Qa,Q/,/i ,h,(Tg,n s ) 



fid 



= (0.27, 0.73, 0.0227, 0.7, 0.8, 0.97) 



(A.l) 



where h = i/rj/100 kms 'Mpc '. These values produce a fidu- 
cial sound horizon of 



Here, we use t he parametrized fitting formula introduced by 
lMcDonaldl d2003) to fit the results of the power spectrum from 
several numerical simulations, 



P F (k,fXk) = b]{\ +^ 2 k ) 2 P L (k)D F (k,fx k ) . 



(A3) 



r sJid = 152.76 Mpc 



(A.2) 



where = k\\/k is the cosine of the angle between k and the 
line of sight, bg is the density bias parameter, /3 is the redshift 
distortion parameter, Pi(k) is the linear matter power spectrum, 
and Df(k,fik) is a non- linear term that approaches unity at small 
k. This form of Pp is the expected one at small k in linear theory, 
and provide s a goo d fit to the 3-D Lya observations reported in 
ISlosar et "aD (f201ll> . We do not generate a density and a velocity 
field, but directly create the Lya forest absorption field instead, 
with the redshift distortions being directly introduced in the in- 
put power spectrum model of equation (IA.3b . with the parameter 
/3 that measures the strength of the redshift distortion. 

To model the evolution of the forest with redshift, b$ 
varies with redshift acc ording to bg - 0.14[(1 + z)/3.25] L9 
dMcDonald et all |2006j). The redshift distortion parameter is 
given a fixed value of 0p _ 1 .4. The non -linear correction factor 
D(k,fXk) is taken from lMcDonald I (120031) . The flux field was con- 
structed by generating Gaussian r andom fields g with an a ppro - 
priately chosen power spectrum (IFont-Ribera et alj 1201 2al) to 
which the log-normal transformation F = exp(-ae vg ) is applied. 
Here a and v are free parameters chosen to r eproduce the flux 
varian ce and mean transmitted flux fraction (McD onald et all 
l2006h . 

DLA's were added to the spectra a ccording to the procedure 
described in lFont-Ribera et all (1201 2bl) . 

Finally, the spectra were modified to include the effects of the 
BOSS spectrograph point spread function (PSF), readout noise, 
photon noise, and flux miscalibration. 

Fifteen independent realizations of the BOSS data were pro- 
duced and analyzed with the same procedures as those for the 
real data. 

We used the mock spectra to understand how our analysis 
procedure modifies the correlation function. Figure I A. 1 1 shows 
the average over 15 mocks of the reconstructed quadrupole and 
monopole using methods 1 and 2 (sections 13.11 and 13.2b and 
that reconstructed with the true continuum. The monopole and 
quadrupole for the two methods have a general shape that fol- 
lows that found with the true continuum including the position 
of the BAO peak. However, both methods produce a monopole 
that becomes negative for 60 fr'Mpc < r < 100 fr'Mpc while 
the true monopole remains positive for all r < 130 Ir'Mpc. As 
discussed in section 13.21 this result is due to the continuum es- 
timation of the two methods which introduced negative correla- 
tions. For both methods, however, the BAO peak remains visible 
with a deviation above the "broadband" correlation function that 
is hardly affected by the distortion. 

Figure IA.2I presents ^o(f) and &(r) found with the data, 
along with the mean of 15 mocks. The figure demonstrates that 
our mocks do not perfectly reproduce the data. In particular, 
for r < 80 fr'Mpc, the monopole is underestimated and the 
quadrupole overestimated. Since we use only peak positions to 
extract cosmological constraints, we only use the mocks quali- 
tatively to search for possible systematic problems in extracting 
the peak position. 

Appendix B: Results for a Fiducial BOSS 
Lya Forest Sample 

The spectra ana lyzed here are all available through SDSS DR9 
dAhnetall l20l2h. and the DR9 quasar catalog is described by 



16 



N.G. Busca et al.: BAO in the Lyo- forest of BOSS quasars 



0.4 



0.2 

0.1 

0.0 

-0.1 
-0.2 



0.3 j% 



Exact Continua 
Method 1 
Method 2 



• * 



V. **** 



50 100 150 

r [h-'Mpc] 



200 



0.4 

0.2 

o 

^ 0.0 

-0.2 
-0.4 



Data (method 2) 
Mean Mocks 



(WW 



50 100 150 

r [h^Mpc] 



200 



0.0 
-0.2 

s ~ 0A 

1h -0.6 

-0.8 
-1.0 



1 .... 1 ... 











50 100 150 

r [h-'Mpc] 



200 



Oh 




-0.5 i 



-1.0 



100 150 
r [h-'Mpc] 



200 



Fig.A.l. The effect of the continuum estimation procedure 
on the correlation function found with the mock spectra. The 
black dots are the average of monopole and quadrupole obtained 
with the 15 sets using the exact continua. The blue (red) dots 
show those obtained with the continuum estimation of method 1 
(method 2) as described in section [3~TI ( 13.2l ). 



Fig. A.2. Comparison of the correlation function for the mock 
spectra and that for the data. The red dots show the mean of the 
15 sets of mock spectra and the black dots show the data. 



iParis et aTl (120121) . A Lya forest analysis requires many detailed 
choices about data selection and continuum determination. To 
aid community analyses and comparison of results from dif- 
ferent groups, Lee et al. (in preparation) will present a fiducial 
BOSS Lye forest sample that uses constrained PCA continuum 
determination (?) and reasonable choices of masks for DLAs, 
BALs, and data reduction artifacts. Both the data selection and 
the continuum determination differ from those used here. Figure 
IB. II compares the Method 2 correlation function from this pa- 
per's analysis to that obtained by applying the Method 2 weights 
and correlation measurement code directly to the continuum- 
normalized spectra of the fiducial Lee et al. sample. The good 
agreement in this figure, together with the good agreement be- 
tween our Method 1 and Method 2 results, demonstrates the ro- 
bustness of the BAO measurement, and the more general corre- 
lation function measurement, to the BOSS DR9 data. 
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r pr'Mpc] 



Fig. B.l. Comparison of the monopole and quadrupole correla- 
tion functions for the sample used here (black dots) and for the 
sample and continua of Lee et al.(in preparation) (red dots). 
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